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ABSTRACT 

This  paper  presents  a  method  for 
real-time  scheduling  and  routing  of  mate¬ 
rial  in  a  Flexible  Manufacturing  System 
(F1!S;  .  It  extends  the  earlier  scheduling 
work  of  Ximemia  and  Gershwin.  The  FMS 
model  includes  machines  that  fail  at  ran¬ 
dom  times  and  stay  down  for  random  lengths 
of  time.  The  new  element  is  the  capabi¬ 
lity  of  different  machines  to  perform  some 
of  the  same  operations.  The  times  that 
different  machines  require  to  perform  the 
same  operation  may  differ.  This  paper 
includes  a  model,  its  analysis,  a  real¬ 
time  algorithm,  and  examples. 

1.  INTRODUCTION 

Purpose 

The  purpose  of  this  paper  is  to  deve¬ 
lop  an  algorithm  to  calculate  real-time 
loading  and  routing  decisions  for  a  Flex¬ 
ible  Manufacturing  System  (FMS) .  An  algo¬ 
rithm  for  calculating  loading  decisions 
for  such  systems  has  been  described  in 
earlier  papers  (Kimemia  and  Gershwin, 

1983;  Gershwin,  Akella,  and  Choong,  1985; 
Akella,  Gershwin,  and  Choong,  1985).  An 
algorithm  for  routing  decisions  is  des¬ 
cribed  in  Maimon  and  Choong  (1985) .  Here, 
routing  and  loading  are  calculated  toge¬ 
ther. 

As  in  the  earlier  papers,  the  problem 
is  to  decide  which  part  should  be  dis¬ 
patched  next  into  a  set  of  machines. 

These  machines  are  capable  of  performing 
work  on  a  set  of  different  part  types  with 
no  time  lost  for  setting  up.  Decisions 
are  made  in  response  to  disruptions  of  the 
operation  of  the  system  caused  by  machine 
failures,  and  according  to  the  surplus  or 
backlog  for  each  part  type.  Whenever  a 
machine  changes  state  (i.e.,  fails  or  is 
repaired) ,  a  new  schedule  and  a  new  rou¬ 
ting  scheme  is  calculated  via  a  feedback 
lav. 


A  limited  form  of  routing  flexibility 
was  allowed  in  the  earlier  work.  Only 
identical  machines  could  perform  the  same 
operation.  In  that  case,  a  part  could  be 
routed  to  the  first  available  copy.  The 
purpose  of  this  paper  is  to  deal  with  sys¬ 
tems  in  which  machines  are  not  identical, 
but  where  different  machines  may  perform 
some  of  the  same  operations.  Different 
machines  may  therefore  have  overlapping 
capability,  and  different  machines  perfor¬ 
ming  the  same  operation  may  take  different 
lengths  of  time  to  do  it.  The  routing 
problem  is  therefore  to  choose  among  al¬ 
ternate  machines  for  some  or  all  the  ope¬ 
rations.  A  model  capable  of  analyzing 
such  issues  model  is  required  for  the 
study  of  certain  real  systems. 

Kimemia  and  Gershwin  (1983)  proposed 
a  routing  algorithm  to  go  along  with  their 
scheduling  scheme.  However,  while  the 
scheduling  method  was  effective,  the  rou¬ 
ting  method  was  not.  In  particular,  the 
routing  decisions  that  would  have  been 
calculated  by  the  method  suggested  there 
might  not  be  feasible. 

Examples 

This  work  was  motivated  by  two  actual 
Flexible  Manufacturing  Systems,  one  from 
the  electronics  industry  and  one  from  the 
metal  cutting  industry. 

We  examined  a  robotic  system  for  the 
assembly  of  printed  circuit  boards  (PCB) , 
particularly  the  part  of  the  system  where 
oddly-shaped  components  are  inserted  tc 
the  board  (such  as  large  electrolytic 
capacitors,  switches,  and  connector 
strips) .  These  components  cannot  be 
inserted  with  existing  dedicated  automated 
machines  (e.g.,  SIP,  DIP,  VCD),  because  of 
their  variability  and  special  handling  and 
assembly  requirements.  However,  some 
types  of  robots  (e.g..  Adept,  IBM  7575), 
equipped  with  appropriate  fixtures  and  end 
effector  tools,  car.  meet  the  job  require- 
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B«nts  (*.9.,  tolerance  better  than 
0.005"),  and  are  adaptable  (progranusable) 

■o  that  they  can  handle  different  types  of 
odd  components. 

As  a  result,  different  operations 
(insertions  of  odd  component  types]  can  be 
performed  by  different  robots,  but  the 
amount  of  time  required  for  a  given  opera¬ 
tion  depends  on  the  speed  of  the  robot 
that  performs  it.  Also,  each  robot  has 
different  configurations  (e.g.,  tools)  and 
inherent  capabilities  (e.g.,  accuracy  and 
repeatability) ,  which  results  in  different 
subsets  of  operations  that  each  rebet  can 
handle  (with  nonempty  intersections  among 
those  subsets) . 

As  a  consequence,  not  only  does  the 
input  rate  of  part  types  into  the  system 
have  to  be  determined,  but  also  the  deci¬ 
sion  of  where  to  send  each  part  for  each 
operation  (among  the  possible  alternative 
robots)  has  to  be  made. 

Such  systems  are  usually  justified 
economically  only  if  the  production  volxune 
is  quite  high  (e.g.,  hundreds  of  thousands 
of  components  inserted  per  year)  and  the 
variety  is  high.  Because  of  their  flexi¬ 
bility,  they  are  expected  to  meet  demands 
that  vary  in  the  short  term  and  that  re¬ 
quire  high  utiliration.  The  work  presen¬ 
ted  here  aims  to  improve  system  perfor¬ 
mance  (e.g.,  to  lead  to  higher  throughput 
and  reduced  WIP  while  meeting  production 
demands) . 

Another  type  of  manufacturing  system 
is  comprised  of  conventional  and  advanced 
machining  centers.  The  latter  are  capable 
of  performing  different  operations,  with 
varied  capabilities.  For  example,  some 
machining  centers  can  do  drilling  and 
milling  operations  that  otherwise  require 
two  different  conventional  machines.  Also 
there  are  3-  and  5-axis  machining  centers. 
The  latter  can  do  more  operations  than  the 
former  without  changing  the  part  fixtu- 
ring. 

As  in  an  electronic  insertion  system, 
the  scheduling  and  routing  problem  in  a 
system  of  several  machining  centers  is  not 
only  to  decide  on  the  input  flow  rate  of 
each  part  type,  but  also  where  each  opera¬ 
tion  should  be  done  among  alternative 
machines  with  different  capabilities. 

■  Literature  Survey 

In  this  paper  we  present  a  method 
that  considers,  at  the  same  time,  two 
functions  --  short-term  scheduling  and 
routing  —  based  on  a  global  view  of  the 
system.  Many  references  consider  just  one 


of  these  functions.  For  example,  Whitt 
(1986)  presents  a  method  which  can  be  used 
just  for  the  local  routing  decisions. 
Although  his  paper  develops  generic 
queueing  methodology,  we  use  his  results 
to  show  an  example  of  local  routing  consi¬ 
derations. 

By  local  routing  decisions  we  refer 
to  a  situation  by  which  a  customer  (or  a 
part)  has  to  join  one  of  several  queues. 
These  queues  represent,  for  example,  the 
input  buffers  to  workstations.  The  alter¬ 
native  queues  are  those  of  the  alternative 
workstations  that  can  perform  the  next 
operation  on  a  part,  which  has  just  fi¬ 
nished  a  particular  operation. 

Whitt  shows  that  in  some  cases,  the 
system  average  delay  is  not  always  mini¬ 
mized  by  customers  joining  the  queue  that 
minimizes  their  own  individual  expected 
delay.  This  result  suggests  that  deci¬ 
sions  should  be  made  only  when  taking  a 
global  view  of  the  system. 

Routing  is  treated  in  papers  by  Hahne 
(1981),  Tsitsiklis  (1981),  and  Seidmann 
and  Schweitzer  (1984).  Hahne  and  Tsitsik¬ 
lis  deal  with  only  two  choices  and  ma¬ 
chines  whose  randomness  is  due  to  failure 
and  repair.  Seidmann  and  Schweitzer  have 
many  choices,  but  the  randomness  is  due  to 
variations  in  processing  times.  Xn  all 
cases,  the  full  system  is  not  considered. 
Instead,  only  one  decision  point  is  con¬ 
sidered,  and  decisions  are  made  on  a  pure¬ 
ly  local  basis. 

By  contrast,  we  consider  the  whole 
system  and  do  not  treat  local  conditions 
in  detail.  This  suggests  that  a  hierar¬ 
chical  decision  policy,  in  which  both 
kinds  of  decisions  —  local  and  global  — 
are  made  separately,  may  be  appropriate. 
The  local  decisions  should  be  made  in  a 
way  that  is  consistent  with  the  decisions 
made  on  a  global  basis. 

Outline  of  Paper 

Section  2  states  the  problem.  Sec¬ 
tion  3  contains  our  solution,  which  is 
based  on  dynamic  programming.  Section  4 
describes  some  numerical  examples  and 
simulation  results.  Conclusions  and  new 
research  directions  are  discussed  in  Sec¬ 
tion  5. 

2.  PRODLEM  STATEMENT 

Section  1  describes  two  situations  in 
which  short-term  scheduling  and  routing 
decisions  are  required.  In  this  section 
we  represent  such  manufacturing  systems 
with  a  mathenatical  model. 


3 


Th«  input  to  th«  problM  is  the  pro¬ 
duction  rsquirsmsnts  and  process  data  in 
the  form  of  process  plans  and  routing 
sheets.  They  specify  the  operations  that 
each  part  type  has  to  go  through,  together 
with  a  partial  precedence  relation  among 
the  operations.  For  each  operation,  a  set 
of  alternative  machines,  and  the  time  for 
the  operation  at  each  machi.ne,  (and  ma¬ 
chine  reliability)  are  specified. 

We  seek  a  feedback  law  which  deter¬ 
mines  when  each  part  should  be  released 
into  the  system  and  which  route  it  should 
take  when  it  enters.  The  release  time  and 
the  route  may  be  functions  of  the  current 
repair  state  of  each  machine  as  well  as 
the  current  production  level  of  each  mart 
type. 

Model 

The  FKS  csnsists  of  K  work  stations, 
and  work  station  m  consists  of  iden¬ 
tically  configured  machines.  A  family  of 
N  part  types  is  being  produced.  The  pro¬ 
duction  rate  of  part  type  n  at  time  t  is 

“n  (  t  )  . 

Let  d„  be  the  demand  rate  for  type 
n  parts.  This  is  a  rate  that  is  specified 
by  higher  level  decision-makers  in  the 
decision  hierarchy.  He  assume  here  that 
it  is  constant  over  the  time  interval  of 
interest.  The  model  is  unchanged  if  it  is 
deterministic  but  time-varying,  but  the 
computation  is  made  more  difficult.  Re¬ 
quirements  are  often  stated  in  terms  of 
production  required  over  some  specified 
time  interval;  we  convert  this  to  demand 
rates. 

Let  x„(t)  be  the  surplus  (if  posi¬ 
tive)  or  backlog  (if  negative)  of  type  n 
parts  at  time  t.  It  is  the  difference 
between  production  and  demand,  and  is 
given  by 

^  -  u„(t)  -  d„.  (1) 

The  states  of  the  work  stations  are 
given  by  a,(t).  This  is  an  integer 
which  indicates  the  number  of  machines  of 
work  station  m  that  are  operational  at 
time  t.  The  vector  a  is  assumed  to  be 
the  state  of  a  continuous  time  Markov 
process  with  rates  K,  so  that 


Let  be  the  rate  at  which 

work  station  m  performs  operation  k  on 
type  n  parts.  (Since  only  a  few  opera¬ 
tions  among  all  those  that  are  possible 
are  performed  on  each  part  type,  most  of 
these  variables  are  0.)  The  relationship 
between  u„  and  y^  is  given 

by  “ 

u-  ■  T  y*  for  any  k  and  n,  (3) 

"  £  m 

In  this  section,  we  formulate  an 
optimization  problem  whose  solution  is  the 
optimal  set  of  y*  variables  as 
a  function  of  time.  In  Section  3,  we 
describe  a  suboptimal  solution. 

Capacity 

The  rate  of  flow  of  material  into  the 
system  is  limited  by  the  rate  at  which 
machines  can  do  operations.  Each  opera¬ 
tion  takes  a  finite  amount  of  time,  and  no 
machine  can  be  busy  more  that  100%  of  the 
time.  A  fundamental  assumption  is  that 
there  is  no  buffering  inside  the  system. 
This  reduces  the  total  work  in  process, 
but  increases  the  need  for  effective  rou¬ 
ting  and  scheduling. 

Let  be  the  amount  of  time  that  a 
machine  in  work  station  m  requires  to  do 
operation  k  on  a  part  of  type  n.  The  rate 
at  which  machines  of  that  station  have  to 
do  such  operations  has  already  been  de¬ 
fined  as  y^  . 

IIB 

During  a  short  interval  of  length  T, 
the  expected  number  of  operations  per¬ 
formed  by  the  machines  is  yJ^T.  (It  is 
assumed  that  the  interval  is  short  so  that 
no  repairs  or  failures  take  place  during 
it.)  The  total  amount  of  time  that  all  of 
the  machines  of  station  m  are  performing 
operation  k  on  part  type  n  is  y^Tj^T. 

The  expected  total  amount  of  time  that  the 
machines  of  station  m  are  performing  all 
operations  on  all  part  types  is 


The  total  amount  of  time  available  on  all 
the  machines  of  station  m  is  a.T  if 
a,  machines  are  operational.  There¬ 
fore, 


prob  [  flt(t-6t)-b  I  a(t)-a  ]  ■  (2) 

Recall  that  different  work  stations 
may  be  available  for  some  operations,  and 
that  they  perform  them  at  different 
speeds.  Routing  is  the  decision  of  which 
work  station  will  perform  each  operation. 


To  summarize,  the  y  flow  rates  must 
satisfy  the  following  set  of  equations  and 
inequalities: 

y‘  z  0  V  k,  m,  n 
na 


(4) 


*  *•  •vary  machine  m.  (5) 


3.  SOLUTION 


S  VlL  ■  I  for  all  k  * 

and  all  part  types  n,  (6) 

where  is  the  name  of  the  first  ope' 
ration  performed  on  parts  of  type  n.  De> 
note  by  fl(a}  the  set  of  all  y  flow  rates 
that  satisfy  (4)  -  (6). 

Note  that  fl(al  is  a  random  set.  As 
machines  fail  and  are  repaired  the  instan¬ 
taneous  capacity  changes.  The  rates  that 
material  flows  into  the  system  must  change 
as  a (a)  changes,  as  well  as  in  anticipa¬ 
tion  of  these  changes.  _  _ 

Cost  Function 

We  seek  a  policy  that  minimizes  a 
cost  of  the  form 

JCx,,  tto,  0)  - 

E  ^  J  g{x(s))ds  I  x(0)-X5,«(0)-aj  j  (7) 

in  which  T  is  the  short  term  period,  such 
as  an  eight  hour  shift  and  g(>)  is  a 
positive  convex  function.  We  assume  the 
cost  function  does  not  reflect  true  costs, 
but  instead  is  chosen  to  lead  to  desirable 
behavior.  Thus,  the  details  of  g(*) 
are  not  important.  In  Section  3  we  de¬ 
scribe  an  approximation  method  which  uses 
only  certain  features  of  the  cost  func¬ 
tion. 

Dynamic  Programming  Formulation 

The  optimization  problem  can  be  writ¬ 
ten: 

minimize  J(x,,  a,,  0) 

subject  to  dynamics  given  by  (1)  and  (2) 
and  y  e  fl (a)  . 

Comparison  with  Kimemia  and  Gershwin 


Following  the  usual  dynamic  program¬ 
ming  practice,  define 

J(x,  a,  t)  - 


u: 


g{x(s))ds 


x(t)«x 


,«(t)-a  J  . 


(8) 


This  function  satisfies  the  Bellman  equa¬ 
tion  (Bertsekas,  1976),  which  takes  the 
following  form: 


_ °  y;^-d„) 

li  D,  t) } .  (9) 

This  equation  has  the  following  in¬ 
terpretation:  we  seek  a  function  J(x,a,t) 
such  that  the  values  of  y(x,m,t)  e  (](a(t)) 
that  minimize  the  right  hand  side  of  (9) 
cause  that  expression  to  be  zero.  This  is 
a  nonlinear  partial  differential  equation 
which  we  cannot  expect  to  have  an  analytic 
solution.  (However,  in  the  case  of  a 
single  part  type  and  a  single  machine, 
Akella  and  Kumar  (1986)  were  able  to  find 
a  closed  form  solution.) 


Xf  (9)  has  a  solution,  the  optimal 
control  y  satisfies  the  following  linear 
programming  problem.  Note  that  the  cost 
coefficients  are  time-varying  . 


(I  " 

subject  to 
yen(e)  ^ 


(10) 


It  is  important  to  recognize  that 
this  is  a  feedback  control  law  since  J  and 
n  are  functions  of  x  and  ot.  The  solu¬ 


Kimemia  and  Gershwin  (1983)  formu¬ 
lated  an  optimization  problem  in  terms  of 
u  of  equation  (3) .  This  formulation  is 
correct  when  there  are  no  route  choices 
except  among  identical  machines.  However, 
they  assumed  that  they  could  Ignore  (6) 
even  when  route  choice  existed,  and  then 
determine  y  from  u  after  solving  the  prob¬ 
lem.  This  assumption  is  not  correct:  the 
above  formulation  is.  Without  (6) ,  the 
choice  of  routes  achieved  may  not  be  feas¬ 
ible,  and  (3)  would  not  necessarily  hold. 


tion  y  is  therefore  a  function  of  x  and 

a  . 

Note  that  J  is  positive  since  it  is 
the  expected  value  of  the  Integral  of  g,  a 
positive  quantity.  Note  also  that  feed¬ 
back  law  (10)  minimizes 

S  ■  *  ti 

while  a  is  constant.  This  is  because  y 
appears  in  (9)  only  in  the  same  term  in 
which  It  appears  in  (11) .  If  a  remains 
constant  long  enough,  and  there  is  a  y  C 
R(a))  such  that  (11)  is  negative,  then  3 


•v«ntually  reaches  a  minimum.  We  call  th^ 
value  of  X  that  produces  this  minimum  the 
hedging  point  and  write  it  x^.  If 
possible  the  production  rate  should  remain 
at  a  rate  that  keeps  x  at  the  hedging 
point.  A  positive  hedging  point  serves  as 
insurance  for  future  disruptions. 

After  J  reaches  this  minimum,  J  and  x 
are  both  constant.  Therefore,  at  the 
minimum. 


Akella,  Maimon,  and  Gershwin  (1987) 
demonstrate  a  technique  for  calculating  a 
set  of  values  for  A(a) ,  b(a) ,  and  c(a), 
from  a  specified  g,  for  a  model  similar  to 
the  one  presented  here. 

4.  EXAMPLES 

Example  1;  Three-Hachine  System 


I  y^  -  d„  .  0  (12) 

and 

||(x|;,«,t)  -  0  (13) 

If  there  is  no  y  e  fl(a)  that  satisfies 
(12),  then  J  cannot  reach  a  minimum  for 
finite  x.  That  is,  the  production  lags 
behind  the  demand  requirements  and  x(t) 
decreases.  This  is  because  too  many  ma¬ 
chines  are  currently  down  to  allow  produc¬ 
tion  to  equal  demand. 

There  are  reasons  to  believe  that  the 
solution  of  linear  programming  problem 
(10)  provides  a  satisfactory  scheduling 
and  routing  algorithm  even  if  an  approxi¬ 
mate  J  function  is  used.  This  was  the 
simulation  experience  reported  by  Gersh¬ 
win,  Akella,  and  Choong  (1985)  and  Akella, 
Gershwin,  and  Choong  (1985) . 

In  addition,  it  is  likely  that  the 
repair  and  failure  processes  are  not  ac¬ 
tually  exponential,  not  actually  indepen¬ 
dent  of  the  machine  utilizations  (as  as¬ 
sumed  in  Section  2 ) ,  and  do  not  have  the 
exact  K  parameters  that  would  be  used  in 
(9)  if  an  exact  solution  could  be  calcu¬ 
lated.  Also,  the  g  function  docs  not 
necessarily  represent  true  costs,  but 
rather  is  chosen  to  obtain  a  desired  beha¬ 
vior.  For  these  reasons,  it  would  be  a 
mistake  to  work  very  hard  to  get  an  exact 
J. 

Therefore,  a  reasonable  strategy  is 
to  select  a  J  function  that  has  the  cor¬ 
rect  qualitative  properties  and  that  is 
easy  to  calculate  and  work  with.  Such  a 
function  is  positive  and  has  a  minimum  at 
the  hedging  point  (for  every  «  such  that 
the  demand  is  feasible  for  that  a) . 
Gershwin,  Akella,  and  Choong  (1983)  use  a 
quadratic  function, 

J  -  |xTA(a)x  +  b(«)\  +  c(«). 


Consider  a  three-machine  system  that 
makes  two  part  types.  Machine  1  can  do 
operations  only  on  Type  1;  Machine  2  can 
only  work  on  Type  2;  and  Machine  3  can  do 
operations  cn  both.  In  fact.  Machine  3 
can  do  the  same  operations  that  Machines  1 
and  2  can  do.  Thus  Type  l  parts  can  go  to 
Machine  1  or  Machine  3  and  Type  2  parts 
can  go  to  Machine  2  or  Machine  3.  The 
problem  is  to  decide  where  to  send  each  of 
the  parts  and  how  frequently  to  send  them 
into  the  system. 

The  capacity  set  Q(a)  is  given  by: 


TiiV},  *  «)  (14) 

ijli,  *  «j  (15) 

^ay'a  *  isy!,  *  “3  (16) 

K’  ° 

The  production  surplus  and  backlog  dyna¬ 
mics  are: 

-  y;  *  - 1.  (18) 

-  y’a  *  yj,  -  **3 


If  J(x,  a,  t)  is  known,  then  the  optimal 
routing  and  scheduling  policy  y  satisfies 

*  yu)  *  *  yL)  ^^o) 

This  is  a  feedback  control  law  since 
the  constraint  set  is  a  function  of  «  and 
the  partial  derivatives  are  functions  of  x 
and  8.  To  solve  this  linear  programming 
problem,  several  cases  must  be  considered. 
Figure  1  demonstrates  the  various  regions 
of  aj/ax-space  that  have  different 
solutions.  The  regions  are  indicated,  as 
well  as  the  values  of  y*^  that 
ar«  optimal  in  thos«  reqions.  A-.sc 
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cated  is  which  of  the  following  conditions 
that  determine  the  values. 

^  >  0  (Regions  I  and  III) 


>  0  (Regions  I  and  II) 

*  yjj  -  0,  -  0  I 

<  0  (Regions  II,  IV,  V,  and  VI) 


^  <  0  (Regions  III,  IV,  V,  and  VI) 

^  <  0  and  ^  >  0  (Region  II) 

I  ®1 

»  -  7^  (E) 

^  >  0  and  ^  <  0  (Region  III) 

-  yia  -  ^ 

If  both  derivatives  are  negative 
(Regions  IV  and  V)  ,  yj^  and 
y^  are  already  determined. 

The  remaining  variables,  yj^  (i  ■  1,  2), 
minimize 

M-yl,  *  Mfz, 

subject  to  (16) .  The  solution  is 

{4  fe  -  ^  V. 

»  y2  •  and  yj  •  0  (G) 

T-J3  '3 

*»  y|3  -  ^  and  yj  -  0.  (H) 

In  each  of  these  regions,  the  control 
y*  moves  the  state  x„  through  the  dyna¬ 
mics  C(18)  and  (19)].  The  state  moves  to 
a  boundary  and  then  to  another  region. 
However,  there  is  one  exception.  In  both 
Regions  IV  and  V  the  state  moves  toward 
the  common  boundary,  which  is  given  by 

If  we  follow  rules  (G)  and  (H) ,  the 
state  will  move  back  and  forth  across  the 
boundary  in  an  unrealistical  manner.  This 
is  called  chatt*rint.  It  occurs 
because  the  problem  is  sinfular ,  and 
a  remed/  is  suggested  'ey  Ger3hv;in,  A/iella, 


and  Choong  (1985).  k  strategy  is  found 
which,  when  x  reaches  Region  VI,  keeps  x 
in  Region  VI.  That  is,  it  maintains  (I). 
It  does  this  by  determining  y^ 
which  minimizes  (21)  subject  to  (16)  and 


This  is  simplified  by  assuming  that  J 
is  epiadratic: 

J  -  |x-A(a)x  +  b(«)\  +  c(«).  (23) 


^  ■  A..X.  +  A..X.  +  b. 


^  -  A3,X-  +  A^jX.  *  bj. 


ix,  ° 


and  y  is  chosen  so  that 


+  bj) 

3^0 


-  (A„x,  +  A.:X;  +  b,)  -  0.  (27) 

3|3 

Since  this  is  true  for  more  than  just 
one  instant,  its  first  derivative  with 
respect  to  t  is  also  0.  That  is, 

[a..*.  +  A,.x.  +  b,: 


-  +  (a,.x,  +  AijXj  +  b,]  -  0,  (28) 

M3 

or. 


•  »  .^,4  »  4*' 
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From  (C)  and  (D) , 

Aa[  ot/T’j+y'^^-d,  ]+Ajj[  tti/TL+y^-dj  ]+b2] 

*23^ 

-+[  A„r  ayr;. +y;3-d,)+A„(a:/T=3+y2^-d,]*b,]-0.  (30) 

NOW  (16)  (as  an  equality)  and  (30) 
are  two  equations  in  two  unknowns,  yj,  and 
y’  .  The  solution  is 


j+«3;Ay-^A.j] 

A.,/' 

-g|J-<*:;A:2+Ajj)+bi+bj 

~  A- 

A.,  t:, 

^:3  ^23  J 

and 

After  X  arrives  at  Region  VI,  it 
stays  in  Region  VI  if  y^  and  are 
given  by  (C)  and  (D)  and  yj^  and  are 
given  by  (31)  and  (32)  .  Chattering  is 
avoided. 

5.  CONCLUSIONS 

This  paper  presents  an  extension  to 
the  earlier  Kimemia  and  Gershwin  work  to 
add  a  real-time  routing  calculation  to 
real-time  scheduling.  Thus  this  model  can 
be  used  for  many  more  types  of  manufactu¬ 
ring  systems. 

Future  work  will  include  the  develop¬ 
ment  of  local  operational  rules  which 
follow  the  system  routing  decisions  calcu¬ 
lated  here,  and  extensive  simulation  of 
various  types  of  industries  to  further 
demonstrate  the  use  of  this  work. 
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